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Abstract 

The scalar and vectorial self energies obtained through QCD sum rules are introduced in the Quantum 
Hadrodynamics (QHD) equations. This QHD and QCD mixing show us that the effect of the density 
on the coupling constants is very small. 

PACS numbers: 21.65.+f, 24.85. +p 

Keywords: Nuclear Matter, QCD sum rules and Self-energies 



I. INTRODUCTION 

The research for the exact equilibrium properties of 
the nuclear matter is an old problem of nuclear 

physics and it continues to have a great interest nowa- 
days. Quantum Hadrodynamics (QHD) || is one of the 
models with larger success to describe the properties of 
the infinite nuclear matter as well as those of finite nuclei. 
QHD describes the nucleon-nucleon interaction through 
the mesons exchange (tt, a, u, p, etc.). Several calcula- 
tions of nuclear structure using QHD and its extensions 
were made with success in the explanation of experimen- 
tal data [[| . For the purposes of this work only the sim- 
plest QHD model will be considered. 

The basic QHD-model that explains the nuclear mat- 
ter includes the nucleon (ip) coupled with sigma (a) and 
omega (u>) mesons. In spite of the pion be the principal 
component of the nucleon-nucleon interaction, it is not 
included because the nuclear matter is an isotropic sys- 
tem with parity conservation. Analyzing the model we 
see that the real part of the scalar term (c-meson) is typi- 
cally of the order of several hundred MeV attractive while 
the real part of the time component of the vector term (u>- 
meson) is typically of the order of several hundred MeV 
repulsive. However the energies involved in problems of 
nuclear structure are only of a few tens of MeV. That 
order of energy is obtained on the QHD model due to a 
large cancellation among the scalar and vector pieces and 
this process of cancelation can be controlled through an 
appropriate choice of the coupling constants. However, 
these so called QHD constants are quite different from 
those ones given in Bonn potential [gj or empirical data 
||. So it is interesting to discuss the validity of these 
coupling constants. 

The success of the Walecka model, which is based on 
Dirac's equation, is due to the mutual cancellation be- 
tween the large scalar and vector potentials. However in 
some works it is discussed that the composed nature 
of the nucleon can suppress the scalar optical potential. 
However, as it is argued in , it can be shown that due 



to the nucleon being immersed in the nuclear media, such 
suppression does not exist. This result makes possible 
the use of the Dirac phenomenology for composed par- 
ticles. Besides, recently the effective field theory (EFT) 
JTo| |l5[ has verified that the QHD models are consistent 
with the symmetries of the Quantum Chromodynamics 
(QCD), the correct theory of the strong interaction. This 
last fact motivates the mixing between QHD and QCD 
accomplished in this work. 

It is known that QCD is the correct theory of the 
strong interaction. However, there is not a perturbative 
treatment for QCD at energies involved in nuclear matter 
problems, so it is interesting to use the nonperturbative 
method given by QCD sum rules that was introduced by 
Shifman, Vainshtein, and Zakharov in the late 1970's jl6| . 
The method consists in describing a correlation function 
in terms of hadron degree of freedom as well as quarks 
degree of freedom. This last one can be written in an op- 
erator product expansion (OPE) where the perturbative 
part (Wilson coefficients) is separate from the nonpertur- 
bative (condensates). The power of this technique is in 
the fact that the nonperturbative part is the same for all 
problems. 

The objective of this work is to get the results obtained 
from QCD sum rules in the media Q and use them on 
QHD equations to analyze the necessary constants to ob- 
tain the saturation point of the infinite nuclear matter. 



II. THE QHD MODEL 

The model used to describe the nuclear matter is from 
Serot and Walecka || and it includes nucleons (ip) in- 
teracting with a and u> mesons. So that the Lagrangian 
density is given by 

C = tpi^df, + g s a - grfup - M)ip + 
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where g s and g v are the mesons coupling constants, 
G>„ — d^u> v — dyU^ is the strength tensor for the vector 
field and M, m s and m v are the nucleon, er-meson and u>- 
meson masses, respectively. The Green function for the 
nucleon is given by 
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Here £( s ) and £j* n represent the scalar and vector self- 
energies, respectively. 

In the more simple model only the contributions of the 
baryons from the Fermi sea will be considered. This is 
equivalent to disregard the contributions coming from an- 
tibarions (Dirac sea). So in agrement with the definition 
of energy-momentum tensor 
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where VEV is the vacuum expectation value of the Too 
and 7 is the spin-isospin degeneracy. The self-energies 
can be obtained through "tadpole" Feynman diagrams 
and are given by the following relationships 
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The expression (Q) should be solved in a self-consistent 
way. Therefore all the ingredients are available to cal- 
culate the energy density. This result is also known as 
Mean Field Theory (MFT), because it can be obtained 
with a mean field approach for the meson fields. 

A simple analysis shows us that the self-energies E/ s ) 
and E(,„) are proportional to the scalar and vectorial den- 
sities, respectively. Then the form presented by Eq. (j^) 
represents an expansion on powers of these densities. On 
the other hand, in agreement with the EFT concepts, the 
energy density is a functional and has an expansion in 
powers of the densities (scalar, vectorial, etc.) organized 



through Georgi's naive dimensional analysis (NDA). So, 
the energy densities have the same form that the expan- 
sions obtained by EFT ||-|||. Therefore, since QHD 
has a foundation on the symmetries of QCD, the next 
step is to calculate the self-energies through the QCD 
methods. However, as the energy levels treated in the nu- 
clear problems are in a nonperturbative regime for QCD, 
the sum rules is the appropriate method to be used. 



III. QCD SUM RULES 



The QCD sum rules method 16 begins with the fol- 
lowing time-ordered correlation function, defined by 

U aP (q) = i fd 4 x e^ x (0\T[ Va (x)fj p (0)} |0> , (9) 



where |0) is the physical nonperturbative vacuum state 
and rj a {x) is an interpolating field with the quantum 
numbers of a nucleon. As proposed by Ioffe [jL7|, the 
proton field is given by 



r)(x) = e a bc (u^CjfiUi) 757^ 
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where u and d are the quark fields, a, 6, c are the color 
indexes and C is the charge-conjugation matrix. Now 
the correlation function can be written as an operator 
product expansion (OPE) whose nonperturbative part 
(condensates) can be separated from the perturbative one 
(Wilson coefficients). The OPE can be generated starting 
from the following expansion for the quark propagator 

BE 
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where {qq) vac is the quark condensate which can be de- 
termined from the Gell-Mann-Oakes-Renner relation, 



2m q (qq) l 



mlfl (1 + 0(ml )) 



In this relation, = 138 MeV is the pion mass, / ff = 93 
MeV is the pion decay constant and m q = (m„+TOd)/2 ~ 
7 ± 2 MeV is the average of the up and down quark 
masses. 

The same correlation function, Eq (^|), can be written 
in a phenomenological way, which mean that the cor- 
relator has a hadronic description through the nucleon 
Green function, Eq (0). The match of the theoretical side 
(OPE) and the phenomenological one is the essence of the 
QCD sum rules. Using the fundamental state of the nu- 
clear matter as vacuum and the interacting propagator 
in the phenomenological side, it is possible to obtain the 
following sum rules m 



2 



m n = -^p{m) PB , 



J (v-QCD) 



64tt z 
UP 



PB ■ 



(12) 
(13) 



Here M represent the borel mass with value near M 
lGeV 2 EtII, and 
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The expression (|12|) is a generalization of the Ioffe's for- 
mula |o] to finite density. The sigma term is estimate 
in Ref p0| as <jb — 45 ± 10 MeV. Taking the ratios of 
the expressions ( |l2| ) and ( 13 ) in relation to the mass Mn 
(Ioffe's formula), they are obtained the following results: 
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Although a self-consistent relation is not obtained as in 
Eq (Q), these expressions have a more fundamental nature 
based on quark and gluon degree of freedom. At this 
point, the idea is to use the expressions (|l6|) and jl7| ) to 
calculate the energy density of the nuclear matter. 



IV. NUMERICAL RESULTS AND CONCLUSION 

Firstly it should be noted that only the ratios between 
the coupling constants and the mesons mass appear in 
the expressions where the energy density of the nuclear 
matter is calculated. Some values for these ratios are 
presented in the table [| In the first line of table | the 
coupling constants were adjusted so that the MFT might 
reproduce the saturation properties of the nuclear mat- 
ter. With the use of these coupling constants, the result 
obtained with MFT are the following (dotted curve on 
Fig. 1): saturation point on Uf = 1.12 fm _1 with a 
energy density by nucleon given by —15.75 MeV. This 
result reproduces the exact equilibrium properties of the 
nuclear matter given by Refs [Qj2| , which is the expected 
result, once the coupling constants were chosen for such 
adjustment to occur. However, these coupling constants 
are very different of the empirical values H. 
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MFT 


3.029 10~ 4 


2.222 10~ 4 


Bonn Potential 


3.233 10" 4 


4.104 10" 4 


Empirical Values 


3.402 10 -4 


3.516 10" 4 



TABLE I. Coupling constants used in MFT and in Bonn 
potential and the empirical values. 
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FIG. 1. Binding energy by nucleon in nuclear matter as 
function of Fermi momentum. In the continuous line is 
adopted the Bonn coupling constants and the self-energies 
obtained with QCD sum rules. In the dotted and dashed 
lines, the original MFT equations are used together with the 
coupling constants from the MFT and Bonn potential, respec- 
tively. 



The most accepted values for these constants in free 
space are found in the Bonn potential ||, for which the 
nucleon-nucleon interaction is adjusted to describe the 
scattering data. The Bonn coupling constants for the 
sigma and omega mesons are given on the second line of 
table | while the empirical values are presented on the 
third line of that table. It must be noted that, in MFT 
the value of the coupling constant for the scalar meson 
is higher than the respective value for the vector meson. 
But in the Bonn potential and in the empirical data there 
is an inversion in the magnitude of these constants. Using 
the Bonn's values in Eq. (|J) , the result presented by the 
dashed curve on Fig. 1 is obtained. That curve shows us 
that the energy density for MFT does not present a good 
behavior. There is no saturation point, in other words, 
there is no formation of nuclear matter. The same result 
is obtained with the empirical constants. The problem is 
in the fact that the use of the Bonn or empirical coupling 
constants implicates in a small increase of the attraction 
but with a much larger increase in the repulsion. It can 
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be argued that the difference between the QHD coupling 
constants and the Bonn values is due to the fact that 
the latter are obtained in free space while the QHD val- 
ues are obtained in nuclear media, but this argument is 
not valid for empirical constants. However, using the 
constants given by Bonn potential, and the self-energies 
found through QCD sum rules, Eqs. ( [l6| ) and (|l7|), with 
(7b = 49 MeV and m q = 6.8 MeV, the result represented 
by the continuous line on Fig. 1 is obtained. That curve 
has a saturation point kp = 1.42 fm _1 with —15.75 MeV 
for the energy density by nucleon, which is the satura- 
tion point of the nuclear matter. Therefore the QCD 
sum rules allows us to reproduce the properties of nu- 
clear matter with the simplest QHD model using more 
realistic values for the coupling constants. 

Since the empiric constants are not so different from 
the Bonn values when they are compared to the QHD 
values, it can be that the effects of the density do not 
alter the values of these constants significantly. On the 
other hand, the exact determination of the parameters 
<tb and m q should allow a better evaluation of the cou- 
pling constants in the media for the various QHD models 
and consequently to determine with more precision their 
real variation in relation to the values of the constants in 
the vacuum. Furthermore this result is another indica- 
tive that the quark degree of freedom takes an important 
place in nuclear problems. 

Finally, it is known that the QCD sum rules work with 
more accuracy when there is high transferred momen- 
tum. So I hope this mix between QCD sum rules and 
QHD model can be applied, with success, in systems in a 
regime of high density and high temperature as neutron 
stars. Besides, in these calculations the contributions of 
the gluon condensate and four quark condensate were not 
included. But it is known that the contribution of these 
terms for the sum rule is very small. However, if more 
precise answers are required, these terms have to be in- 
cluded as well as the most sophisticated versions of QHD. 
These and other questions are left for future works. 
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